CASCADE OF PHASE SHIFTS AND CREATION OF 

NONLINEAR FOCAL POINTS FOR SUPERCRITICAL 

SEMICLASSICAL HARTREE EQUATION 
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^\ ' Abstract. We consider the semiclassical limit of the Hartree equation 

("^ , with a data causing a focusing at a point. We study the asymptotic 

(«_») ' behavior of phase function associated with the WKB approximation 

Cn I near the caustic when a nonhnearity is supercritical. In this case, it 

Q_) . is known that a phase shift occurs in a neighborhood of focusing time 

O ' in the case of focusing cubic nonlinear Schrodinger equation. Thanks 

[] , to the smoothness of the nonlocal nonlinearities, we justify the WKB- 

^^ ■ type approximation of the solution for a data which is larger than in 

^^ ' the previous results and is not necessarily well-prepared. We also show 

I |, by an analysis of the limit hydrodynamical equaiton that, however, this 

PLh ' WKB-type approximation breaks down before reaching the focal point: 

.^^ I Nonlinear effects lead to the formation of singularity of the leading term 

• . of the phase function. 

^- 
-(— > . 

C^ . 

S: 

'""'■ 1. Introduction 

\^J I This paper is devoted to the study of the semiclassical limit e ^ for 

,_^ ■ the Cauchy problem of the semiclassical nonlinear Schrodinger equation for 

(N '. {t, x)eR+x W 

^- 2 

^. ■ (1.1) iedtu'' + —Au'=e''N{u'), uIq{x) =ao{x)e-'^-^ 

^^ . with the nonlocal nonhnearity of Hartree type 

O: (1.2) N{u') = X{\x\-^*\u'fX, 

.!_^ I where n ^ 3, A is a real number, and 7 is a positive number. In this paper 

^ ■ we consider the small-nonlinearity case a > 0. The case a = is studied 

in [1, 9]. In the case of linear equation A^ = 0, the quadratic oscillation in 
the initial data causes a caustic at the origin at f = 1. In [4, 5], R. Carles 
justified the general heuristics presented in [20] in the case of (1.1) with 
^{y) = \y\ y (see, also [7]), and two different notions of the criticality for 
a are realized. One is concerned with the nonlinear effect on the behavior 
far from the focal point, and the other with that near the focal point. The 
situation is similar in the case of Hartree equation 

(1.3) iedtu" + —Au' = Ae°(|a;r^ * |n^P)n^ 

(1.4) W|i=o(^)=ao(x)e-*^. 

(see, [8, 10, 27]). The two critical indices are a = 1 (far from focal point) 
and a = 7 (near the focal point). They are completely different notions. For 
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example, the second index = 7 depends on the shape of the nonhnearity 
while the first not. Using terminologies in [6], we have the following nine 
cases: 





a > 1 


a = 1 


a < 1 


a > 7 


linear WKB 
linear caustic 


nonhnear WKB 
linear caustic 


supercritical WKB 
linear caustic 


a = 7 


linear WKB 
nonlinear caustic 


nonhnear WKB 
nonlinear caustic 


supercritical WKB 
nonhnear caustic 


a < 7 


linear WKB 
supercritical caustic 


nonlinear WKB 
supercritical caustic 


supercritical WKB 
supercritical caustic 



Roughly speaking, the term "linear WKB" means that, far from the caustic, 
the propagation of u^ does not involve the nonlinear effect at leading order. 
The term "linear caustic" means that the nonlinear effect is negligible at 
leading order when the solution crosses the focal point. 

Now, let us be more precise about this problem with (1.3)-(1.4). Let w^ 
be the solution of the linear equation {iedt + (e^/2)A)ti;^ = with the same 



initial data as in (1.4). 
approximate solution 



Note that w'^ is approximated by the WKB-type 



(1.5) 



<in 



(1 - t)«/2 



Co 



l-t 



2e(t-l) 



as long as e/(l — t) is small. In the linear WKB case a > 1, it is shown in 
[8, 26, 27] that the solution u^ is close to w^ (and so to ffjjj before caustic. 
The time interval in which n^ is approximated by w"^ depends on the latter 
critical notion = 7: If a > 7 then u^ ^ it;^ as e ^ for all t ^ 1; if 
a = 7 then it holds for 1 — t ^ Ae with a large A; if a < 7 then it holds for 
1 — t ^ Ae^ with a large A and some ^ = ^(a, 7) ^ (a — l)/(7 — 1). In the 
case a ^ 7, the asymptotic profile of the solution beyond the caustic is also 
given in [8]. Let us proceed to the supercritical caustic case a < 7. In this 

case, infinitely many phase shifts occur between the time 1 — t = AieT-i , 
called a first boundary layer, and the time 1 — t = Ae"'''', called a final 
layer. This phenomena, called cascade of phase shifts, is first shown in [6] 
for (1.1) with certain class of nonlinearities including the cubic nonhnearity 
^{y) = \y\^yi ^-^d the asymptotic behavior of the solution is given before 
the final layer, 1 — t S> £"^1/7-1^ jj^ the case (7 >)a > 1. The similar result 
is proven also in the nonlinear and supercritical WKB case a ^ 1, provided 
the initial data is a properly-modified one of the form 



(1.6) 



'■|t=ol 



60 (e' 



,x e 



' 2e 



exp(JeT (l)Q{ei ,x)), 



where (6o(i; x), (j)o{t, x)) is a suitable function defined in terms of oq- Let us 
call this type of data as a well-prepared data. 

The aim of this paper is to give an explicit asymptotic profile of the 
solution of (1.3)-(1.4) (before and) on the final layer for whole 7 > a > 
with a not-modified data. In particular, the behavior of the solution is new 
in the following situations: 

• On the final layer, that is, at t = 1 — T^^e"''^ with not necessarily 
small T > 0. 

• a ^ 1 with an initial data which is not necessarily of the form (1.6). 
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Moreover, it will be shown that the WKB-type approximation of the solution 
breaks down at some time on the final layer, that is, at t = tc := 1 — 
{T*)~^e'^''^ for some T* > for a certain class of initial data. 

For our analysis, we apply the following transform introduced in [6]: 

1 I £ ^^ X \ ( \x\ 

(1.7) u=(t,a;) = — — -— ^M— — ,— — exp i 



(l_t)n/2-^ \X_i^X_tj ^V2e(t-1) 
Then, (1.3)^(1.4) becomes 

1 — — o 

^e^-%^^^^ + II—Ia.V'^ = \T^-Wv\-^ * \r\'W, i^ir=s^/.iy) = «o(y), 

where r = e°''^ /{I — t) and y = x/{l — t). The quadratic oscillation of the 
initial data is canceled out. This transformation clarifies the problem; if the 
nonlinear effect is weak and so if u^ behaves like w^ and vf-^^, then ip^ is close 
to oq; if the solution has a rapid oscillation other than exp(i|xp/2e(t — 1)), 
then ^p'^ becomes oscillatory. We change the parameter into h = e^""/'^. 
The limit e ^ is equivalent to /i ^ as long as a < 7. Denoting ijj^ by 
ip^, our problem is reduced to the limit /i ^ of the solution to 

(1.8) i/^5,V'" + ^A,V" = ArT-2(|y|-'^*|V^'^|2)^^ # ^^(v) = ao{y). 

Since r = h°''^''~°'' /{I — t), the correspondence between boundary layers of 
t and r variables is as follows: 

initial time: t = < > r = /if-" , 

(1-9) first layer: t = 1 - Aie^ < — >t = A^^/i^V^, 

—1 — 
final layer: t = 1 — T sf < — > t = T. 

Our analysis of (1.8) is based on a generalized WKB method by Gerard 
[15] and Grenier [18]. We apply a modified Madelung transformation 

(1.10) ip'' = a'' exp (i^ 

to (1.8), where a^ is complex valued and <p^ is real valued. This choice is 
slightly different from the usual Madelung transformation 

h / h I ■ 

ip = V /^ 6Xp I-— 

V h 

(see [14]) which leads us to an equation of compressible fluid with the quan- 
tum pressure. It is essential that a^ takes complex value and, therefore, 
^h _/. gh jj^ general. The choice (1.10) allows us to rewrite (1.8) as the 
system for the pair (a'*, (j)^): 

dra^ + V0^ • Va^ + -a'^Ac/)^ = i-Aa'^, 



(1.11) 



0. 



a!" c = 


= oo, 


<i)^ 


r=/iT-Q 




\T=ht-'^ 
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If we approximate the solution of (1.11) up to h^ order, that is, if we estabhsh 
an asymptotics such as 

(1.12) a'' = bo + hbi + o{h^), (t)^ = (t)Q + h(t)i+o{h^), 

then, by means of (1.10), we immediately obtain the WKB-type approximate 
solution of if)^. This method is first employed for nonlinear Schrodinger 
equation with a certain class of defocussing local nonlinearity including the 
cubic nonlinearity N{y) = \y\^y for analytic data [15] and for Sobolev data 
[18] (see, also [2, 11]), and is extended to the Gross-Pitaevskii equation [3] 
(see, also [22]) and to the equation with (de)focusing nonlocal nonlinearities: 
Schrodinger-Poisson equation [1, 23] (see, also [21, 30]) and Hartree equation 
[9]. 

Our goal is to justify (1.12). The natural choice of the main part {ho, (po) 



may be 

system 



(1.13) 



^h Mh\ 



\h=Q- 



Letting h = in (1.11), we obtain a hydrodynamical 



drbo 


+ V(/.o- 


V6o + -h 


3A(/.o 


5r0O+2lW 


ol' + Ar^- 


-\\y\ 


^0|r= 


=0 = ao, 


0O|r=O = 


0. 



-7 



0, 



For a Sobolev data oq, (1.13) has a unique local solution (Theorem 1.3). 
Then, the main task is to determine /i^-term (6i,(/>i). The difficulty of 
finding /i^-term lies in the following two respects: 

(1) The equaiton (1.11) itself depends on h through the term i-^Aa^. 

a. 

(2) The initial time of (1.11) tends to r = at a speed h'<-" . 

We will see that the first becomes crucial when we consider the asymptotic 

behavior of the solution on the final layer, and that a suitable choice of 

h^-term is the key for overcoming this difficulty. We use (fteqm ';^equ) defined 

by (1.18), below, as an /i^-term and show an asymptotic behavior of t/j^ for 

T G [hi-" ,T] when a ^ 1, where T is independent of h (Theorem 1.3). By 

(1.9), T G [/i^,r] is equivalent to t G [0, 1 - T-^e^/T], from the initial 

time to the final layer. On the other hand, the second becomes crucial in 

the supercritical case a < 1. This is because the moving speed of initial time 

becomes too slow. In this case, we need three more kinds of correction terms 

whose order are between /i" and h^. With them, we describe the asymptotic 

behavior also for a g]0, 1[. A heuristic observation on these correction terms 

is in Section 2.2. The rigorous result is in Theorem 5.1. 

The well-prepared data (1.6) is closely related to the second problem listed 

above. The function (6o,</'o) in (1-6) is the solution of (1.13). If we employ 

the well-prepared data and consider (1.3) with (1.6), then (1.11) changes 

into 

I h . ,h .h 

—( 

2 



(1.14) 



dra^ + V(P^ • Va^ + 4a^ 



^-Aa^' 



1 



T=h-f-°' T = 



* a' 



h\2\ 



0, 



_c^ = (t)o{h~ 
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The initial time is still moving, however, the only difference between (1.14) 
and (1.13) is the existence of ii^^a ■ Thus, we will see that we do not meet 
with the second problem any longer. This point is discussed in Sections 5.2. 
We also consider the problem of global existence of the solution to (1.13). 
With the notation {p,v) := (|6oP) V(/)o), (1.13) is the compressible Euler 
equation with time-dependent pressure term of Hartree type: 



(1.15) 



' drP + div(pt;) = 0, 
drv + {v ■ V)f + ArT~2v(|xr^ * p) 

,/>|r=0 = |aOp, W|r=0 = 0. 



We adapt results in [12, 25] (see also [24, 28]) to prove that C^-solution to 
(1.15) cannot be global in time, in several situations (Theorems 1.5 and 1.6). 
This implies (1-12) breaks down before caustic t = 1 (see Section 2.1). 

1.1. Main result I. To state our result precisely, we introduce some nota- 
tion. For n ^ 3, s > n/2+ 1, p G [1, oo], and q G [1, oo], we define a function 
space yp%(M") by 



oo/ 



(1.16) y;^^{w^) = c- 

with norm 

(^ ^7\ 11.11 •— 11.11 4- ll\7.ll 4- ll\7^.ll 

\^-^' ) II lly/_5(M") •" II IILP(R") ^ II ^ llL9(Rn) T^ II ^ llH«-2(Rn) " 

We denote llf = K* (M"), for short. For q < n, we use the notation 
q* = nq/{n — q). This space Ypn is a modification of the Zhidkov space X'^, 
which is defined, for s > n/2, by X'{W') := {/ G L°°(M")|V/ G H^~^{W)}. 
The Zhidkov space was introduced in [31] (see, also [13]). Roughly speaking, 
the exponents p and q in Y^ indicate the decay rates at the spatial infinity 
of a function and of its first derivative, respectively. Moreover, the Zhidkov 
space X^ corresponds to Y^ 2 '^^ ^ sense, if n ^ 3. We discuss these points 
more precisely in Section 3.1. We also note that Y22 is the usual Sobolev 
space H^. We use the following notation: Y^ := r\s>oY^^,q] for intervals Ii 
and I2 of [l,cx)], Yj\^^ := Hpei^Y^^g and Y^j^ := Hge/j^p'g. These notation 
are sometimes used simultaneously, for example Y^j^ := ^s>o,peh,qei2^p,q- 
We also use the operator 

^,2 . |2 

This is the scaled version of the Galilean operator and is suitable for a study 
of rapid phases other than e*'^' /2e(*-i)_ 

We also introduce the following systems for a pair {bcqu,4>cqu)'- 
(1.18) 

111 

drbcqu + V0cqu " V6o + VcpQ ' Vfecqu + -bcqu^4'0 + -^oA^cqu = 9^^°' 
^drCpcqn + V(/>0 • V0cqu + At'^"^ (|y | "'^ *2Rebobcqu) = 0, 

where (bo, (po) is a solution of (1.13). This will be posed with the zero data 

(1-19) &equ|T=0 = 0> <Aequ|T=0 = 
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or the data 

(1.20) fccqu|r=o = 0, 0cqu|r=o = A(|yr'^ * |ao|^)/(7 - !)• 

Notation 1.1. Let T > and X be a Banach space. Let {%} be an increasing 
sequence of real number, 4>{t,x) £ C{[0,T];X) be a function , and {(pj} be 
a sequence of function in X. We write 

oo 

(p{t,x)>i'^t''^(l)j inX 



if it holds that 



as t ^ for all J > 1. 



0(t,x)-^t'=^ 



i=i 



ait'-'] 



X 



We now state our main result. To avoid complicity, here we state the 
result for a ^ 1; the linear and nonlinear WKB cases. For the supercritical 
WKB case a < 1, see Theorem 5.1. 

Assumption 1.2. Let n > A and A G M. The constants 7 and a satisfy 
max(l,n/2 — 2)<7<n — 2 and < a < 7, respectively. The initial data 
ao G H°°. 

Theorem 1.3. Let assumption 1.2 he satisfied. Assume a ^ 1. Then, there 
exists an existence time T > independent of e. There also exist (6o,<Ao); 
(6equ,<^cqu) G Ci[0,T];H^ X l^(^/,,^],(„/(,+i),^] ) such that: 

(1) Mr,y) ^ Ej°^i^^^"Vi(y) in ^(^/^,oo],{n/(7+i),oor 

(2) The solution u^ to (1.3) with (1.4) satisfies the following asymptotics 
for all s ^ 0: 
(1.21) 



sup 

te[0,l-T-ie"/T] 



|J"|Mn^(t)e-^*'W 

as e ^ with 

(1.22) $=(t,x) =e^~^(^ 
and 

(1.23) A'{t,x) = bo 



2 



(1 - t)«/2 ^ ' 







eT 



x 



1-t' 1-t 



eT 



ET 



where {bQ,WQ) solves (1.13) and (ftcqu, ';^cqu) solves (1.18) with (1.19) 
if a > 1 and with (1.20) if a = 1. 

Remark 1.4. (1) By the definition of "x" sign, the expansion (po{T, y) x 

Ei^iT'^^~Vj(y) implies 
J 
(1.24) Mr) = E^'^'V, + o{r^''') in l(^/,,oo],(„/(,+i),^] 

i=i 
as r ^ for all J > 1. 
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(2) In Theorem 1.3, we only need uq G H^" with some sq > n/2 + 
3. Then, {bo, <Po) belongs to C([0, T],H^o ^ y^^+^^^ ^^/^^^^ ^,) and 

(6equ,<Aequ) bdongS tO C{[0,T], H^O-2 ^ Y^^;^/^^^, („/(^^,)^^,). The 

asymptotics (1-21) holds for any s £ [0, sq — A]. 

(3) All ifj are given explicitly (but inductively) in terms of oq. For 
example, ifi = A(|x|~''' * |aoP)/(l — 7) and 

^^=- 2(7-1^27-1) '^^'"''"*'"°''^'' 

- ,(^ _ i)(2^ _ 1) i\-r * (V • (|aopV(|.|- * Worn 

(see. Proposition 4.2). 

(4) Even if the system (1.18) is posed with zero initial condition, its 
solution (6equ, i;^equ) IS not identically zero because of the presence of 
the (nontrivial) external force |A6o- 

(5) The choice of the initial data (1.20) is the key for the analysis in the 
case a = 1. We discuss this point more precisely in Section 2.2. 

The asymptotics (1.21) reads 

(1.25) u%t, x) ~ ^-±-^A%t, x)e^*^(*'-)e^^ 

as e ^ 0. Indeed, this holds in -L°°([0, 1 - r-^e^/T]; L^). This explains the 
cascade of phase shifts. We consider the case 1 < a(< 7). Combining (1.24) 
and (1.22), we have 

■^ gctj-i / X \ / e""'"^ \ 

We set g^Jt,x) = n-fV-v-i '^j(l^)- Then, the above asymptotics yields 



{i-ty 
$=(t,x)~0 forl-t>e^F^, 



a-l 



2c«-l 



<^%t,x) r^ gl{t,x) for 1 -t > £2^-1, 

<I>^(t,x) r^ gl{t,x) +gl{t,x) for 1 - t >e3f^. 



Jg-l 



J 

<I>^(t,x) ~ ^g|(t,x) forl-t>e7^^ 



as e ^ 0. On the other hand, the amplitude A^ satisfies 

A^{t,x) = aQ(-^\ +0(1) 



8 S. MASAKI 

as ST /(I — t) — > 0. Substitute these expansions to (1-25) to obtain 

u^ ~ vL = (TTTj^^o (^47) ^'^^ for 1 - * » ^^' 

n^ ~ t;fi„e^3i(*'^) for 1 - t > e^ , 

u' ~ ^;f.^eis!(t:^)+isi(t,^) for 1 - t > e^^ , 



u^ ~ ufi^e* ^J=i ^l*^*'""^ for 1 - t > e J7-1 , 



Recall that vf-^^, given in (1.5), is the approximate solution for the linear 
solution w^. One sees that the solution behaves like a free solution in the 

g-l 

region 1 — t ^ st-i where the initial time t = lies, and that, at each 

boundary layer of size 1 — t ~ e -^t^i (the J-th boundary layer), a new phase 
associated with Qj becomes relevant. 

1.2. Main result II. Our next result is the non-existence of a global so- 
lution to (1.13). We further assume the radial symmetry and 7 = n — 2 
(n ^ 3) in (1.15). A suitable change of A yields the radial compressible 
Euler-Poisson equations 

( dripr""-^) + dripyr""-^) = 0, 



(1.26) 



drv + vdrv - Xr^'-^drVp = 0, 

a,(r"-i5,yp) = pr"-\ 

, P\t=o = koP, 1'|t=0 = 0, 



where r := |x|. We define the "mean mass" Mq in {|x| ^ r} by 
Mo{r):=^ r \ao{s)fs^-'ds. 



r" JO 



Theorem 1.5. Let A < and n ^ 4. For every nonzero initial amplitude 
Oo G C^ , the solution to (1.26) breaks down no latter than 



\>^\sup.,^QMo{r) 

Theorem 1.6 ([12], Theorem 5.10). Let X > and n = A. The radial C^- 
solution (a, f) to (1.26) is global if and only if the initial amplitude oq G C^ 
satisfies 

|ao(r)|2 ^ 2Mo(r) = 4 / \aois)\'^ s^ ds 
r Jo 

for all r ^ 0. In particular, if aQ G L^(R*^) then the solution breaks down in 

finite time. Moreover, the critical time is given by 

1 

2 



Amaxr>0 (2Mo(r) - |ao(r)|^ 
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Example 1.7 (An example of finite-time breakdown). Consider the equation 
(1.26). Let n = 4 and A > 0. Suppose 

ao{x) = ao{r) = r '^e a^ . 

Then, an elementary calculation shows 

Mo(r) = r'^^e'^ 

and so that 2Mo(r) — |aoP = (2?^ — l)r~^e^~ takes its maximum at tq = 
(7 + ^/l7)/16 = 0.6951 . . . which is the root of 

Sri - 7ro + 1 = 0. 

By Theorem 1.6, the critical time is 



5 / 3405 + 827^/17 3(7-vT7) 
— e 32 



JA|(2Mo(ro)-|ao(ro)|2); V A2i3 

In fact, the solution to (1.26) is given by 

\a\\T,X{t,R)) ^^' 



X2(t, R){2R^eR - \{2R - 1)t^) 
v{T,X{r,R))= ^l 



where 



X{t, R) = Ryi + XR-^e-RT^. 



At the time t = Tc, the characteristic curves "touch" at r = re := X(ic,ro), 
that is, we have (9i?X)(Tc,ro) = 0, which is one of the sufficient and nec- 
essary condition for finite-time breakdown (see, [12]). More explicitly, we 
can see that, as r tends to Tc, the amplitude |ap blows up at re since the 
denominator 

2r^e^ - A(2ro - l)t^ 

tends to zero as r ^ Tc. We illustrate the calculation in Remark 6.2. 

In this example, oq G H°° C C^ and so Theorem 1.3 holds. We see that, 
however, the asymptotics (1.21) is valid only for s'< /{I — t) < Tc and it 
cannot hold for e^ /(I — t) ^ Tc- 

The rest of the paper is organized as follows. In Section 2, we make a 
summary of the results in this paper with previous results. Section 3 is 
devoted to preliminary results. We prove Theorem 1.3 in Section 4. The 
strategy of the proof is illustrated rather precisely in Section 4.1. In Section 
5, we treat the supercritical WKB case a < 1 (Theorem 5.1). The well- 
prepared data is discussed in Section 5.2 Finally, we prove Theorems 1.5 
and 1.6 in Section 6. 
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2. Summary 

2.1. Cascade of phase shifts in the hnear WKB case. We first discuss 
about the linear WKB case a > 1. According to Theorems 1.3, 1.5, and 
1.6, we summarize the result in this case as follows. Recall several boundary 
layers: 

initial time: t = < > 1 — t = 1, 



(2.1) 



a—1 a—1 

first layer: t = 1 — Aist-i < — > 1 — t = Aict-i , 
J-th layer: t = 1 - Aje -^-'-^ < — > I -t = Aje ^^-i , 



final layer: t = l-T"^e^ < >l-t = T~^£-', 

where Aj and T are positive constants. 

• From the initial time to the first layer. Before the first layer, 

that is, for 1 — t S> e-'-'^ , the behavior of the solution is the same as in the 
linear case at leading order. Indeed, the asymptotics (1.21) implies that the 

solution behaves like vf^^ defined in (1.5) for 1 — t S> e''^^ . The phase shifts 



disappear: Prom (1.22) and the time expansion 0o(t) ^ J27Li'^^'' ^V'ji ^^ 



see 

oo / 2L \ 7i-l 

e-i 



|$'(t)|li^(Kn) ^C-e^ i^f__j ||99j-||^^(^„^<l 



a-l 



t 



if 1 — t ^ gT-i . Moreover, the amplitude tends to a rescaling of oq: By 
(1.23), 

A^t, x)^bo( 0, YZ^ j exp i i<pcqn i 0, j^^ j j ="0(7 

since e°''^ /{I — t) <C 1 for 1 — t ^ e^-i . This agrees with the analysis in 
[26]. On the other hand, on the first layer 1 — t = AiEt-i the nonlinear 
effect becomes relevant. The term T'^~^ipi in the sum-*^ Xl^i t"''^~^^j is no 
longer negligible. 

• From the first layer to the J-th layer. Soon after the first layer, 
the solution becomes strongly oscillatory by T^^^tpi. Between the first and 
the second layers, only r^^^tpi is effective because 

/ a \ 7-1 



i-t ' \i-t 



L°°(R") 



7i-i 



— — 1 V"^ IE"' \ 

2a-l 2a-l 

holds if 1 — t S> e ^''^^ . When we reached to the second layer l — t = A2e ^7-1 ^ 
the phase t^'''~^c/?2 become relevant. Similarly, at each J-th layer 1 — t = 

Jg-l 

Aje-'T-i new phase (fj becomes relevant. This is the cascade of phase shift. 



This sum is the formal one: Here, we say "sum" in the sense that 4>o{t) 
J=l ' 



EOO "/ 7 — 1 
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In this regime, e^/(l — t) converges to zero as e ^ 0, and so that the 
ampUtude stih stays the hnear one. 

• On the final layer. After a countable number of boundary layers, we 
reach to the final layer. In this layer, ei /(I — t) does not tend to zero any 
longer. Therefore, the asymptotics of amplitude changes into a nonlinear 
one. It turns out that the ratio T := e°'''^ /{I — t) plays the crucial role. If 
T ^ is small, then the asymptotic behavior of the solution is described 
by the asymptotics (1.21). It is essential to use <I>^ and A^ defined in (1.22) 
and (1.23), respectively. Thanks to the nontrivial remainder term given by 
(/>equ, (1-21) gives the asymptotic behavior of the solution on the final layer. 

On the other hand, (1-21) breaks down at T = T* < cx) in several cases 
(Theorems 1.5 and 1.6). It is because the nonlinear effects cause a formation 
of singularity. These focal points are moving: If the focal point is 

/ — \ 



yi-tci-tc^ 

then, in the (t, x)-coordinates, this focal point is 

which tends to (t,x) = (1,0) as e ^ 0. Example 1.7 gives an example of 
this type blowup. 

• After the final layer. It remains open what happens after the break- 
down of (1.21). However, we can at least expect that more rapid oscilla- 
tions might not appear after the final layer. This is because, in the case 
of A > 0, 2 ^ 7 < 4, and a g] max(7/2, 1),7[, the order of the upper 
bound of ||J^u^||^2 stays e^~"/T even after the caustic ([26]). Recall that 
j£ — ^i\x\ /2£(t~i)^f^^_-^^y^-t\x\ /2£(t-i) gj|^g],g Q^^ ^]^g main quadratic phase. 

Therefore, this divergent upper bound implies that the order of magnitude 
of energy of the oscillation other than the main quadratic phase is at most 

g.l-a/7_ 

2.2. Cascade of phase shifts in the nonlinear and supercritical WKB 
case. We now turn to the case a < 1, the supercritical WKB case. Our 
goal is to obtain a WKB-type approximation similar to (1.21), which ex- 
plains the cascade of phase shifts and gives the asymptotic profile of the 
solution before and on the final layer. We do not use modified the initial 
data (1.6) and keep working with (1.3)-(1.4). As stated above, the analysis 
of (1.3)-(1.4) is reduced to the analysis of (1.11) via the transform (1.7). 
Recall that h = e^~"''^. The difficulty is the following: 

(1) The equaiton (1.11) itself depends on h through the term i-^Aa^. 

(2) The initial time of (1.11) tends to r = at a speed h''-" . 

The second is the main point of the supercritical case a < 1. In this section, 
we discuss by heuristic arguments what happens, what is the problem, and 
how we can overcome it. Then, it turns out that the situation becomes 
more complicated, and the cascade of phase shifts phenomena involves more 
phase shifts and boundary layers than the linear WKB case. The rigorous 
result is in Theorem 5.1 (see, also Remark 5.3). 
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• The transpose of the initial time and the first boundary layer. 

We see from (2.1) that the relation a < 1 causes the transpose of the initial 
time t = and the first boundary layer, that is, the initial time t = 

a— 1 a— 1 

lies beyond the first boundary layer t = 1 — Aict-i because Ai^t-i ^ 1 
for small e. It means that there is no linear regime and the behavior of 
the solution involves nonlinear effects at leading order soon after the initial 
time. 

• The nonlinear behavior of the solution at the initial time. As 

a matter of fact, the behavior of u^ is already "nonlinear" at t = 0. More 
precisely, the principal part of the phase shift of the solution u^ other than 
exp(i|a:p/2e(t — 1)) (we call this as a principal nonlinear phase of u^) for 
a < 1 is given by 

/ / - 

(2.2) exp ieT cf) 




1-t'l 



where (po is the solution of (1.13). Here we remark that the principal non- 
linear phase is the same as in the case a > 1 (In Theorem 1.3, we denote 
this by exp(i<I>'^)). Since (po is given as a solution of (1.13), the shape of (po is 
completely independent of a. Moreover, the choice of (/)q is natural because 
the function (pQ is the unique limit of the phase function </> which solves 
(1.11) with a . Recall that ip = a exp(i(/) //i) is an exact solution of (1.8). 
Using the expansion (Po{t) x J27LiT'^''^^y^j{y)j we have 

/ 2L \ / / « \ 7-1^ 



as long as e~/(l — t) <^ 1. If the right hand side is small, then the phase 
shift caused by nonlinear effects is negligible and so the solution behaves like 
the linear solution. However, the right hand side is 0(e"~^) at t = 0, which 
is not small if a < 1. In this sense, the behavior of u^ is nonlinear at t = 0. 

• The initial condition as a constraint. On the other hand, we always 

have M^(0,x) = ao(x) exp(— i^) since, as stated above, we do not modify 
the initial data and keep working with (1.4). This initial condition seems 
to be quite natural and the simplest one for this problem. This is true if 
a > 1. However, in the supercritical case q < 1, the meaning of this condi- 
tion slightly changes and this condition becomes a sort of constraint: The 
appearing nonlinear effects must disappear at i = 0. To achieve this con- 
straint, we need to employ some more nontrivial phase shifts as correction 
terms in order to cancel out the nonlinear effect at t = 0. This modification 
with correction terms is the heart of the matter. The main difference be- 
tween two initial data (1.4) and (1.6) is this point. A use of (1.6) enables 
us to leave the behavior of the solution at the initial time nonlinear. Hence, 
we do not need any correction term. 

• A formal construction of correction terms. Intuitively, this mod- 
ification is done as follows: First, we replace (/)o(r, y) in (2.2) by (po{T,y) — 
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^ • j^ T^^ ^ipj{y) for some k ^ 2. This yields the following modified princi- 
pal nonlinear phase 



(2.3) exp ie~ 




At t = 0, it holds that 






ffT 



fc-i ._^ 



0{e 



ak-l\ 



Therefore, if we take k large enough then the modified approximate solu- 
tion which the above principal phase (2.3) gives possesses the desired two 
properties: The leading term of the principal nonlinear phase is the same 
as (2.2); and the phase shifts tend to zero at the initial time t = 0. Of 
course, this simple modification is too rude and so valid only in the small 
neighborhood of t = 0. Hence, to obtain an approximation also outside the 
small neighborhood, we replace each ipj{y) by a function —(j)pha.,j{T, y) which 
solves a kind of j'-th linearized system of (1-11) with — </'pha,j(0) = <Pj- This 
yields the principal nonlinear phase like 
(2.4) 



exp zeT 



00 



ST 



l-t l-t 



+ 



E 



c 



yj- 



'f'phaj' 



teT 



l-t' l-t 



We note that teT/(l — t) = ei- /(I — t) — e'l , which is zero at t = 0. 

• Three kinds of correction terms. In fact, the above modified prin- 
cipal nonlinear phase (2.4) is still insufficient. We need two more kinds of 
correction terms which are essentially different from (j)pha,j- Now, let us list 
all kinds of correction terms which we use: 

(1) Correction from phase. The first one is the above (pphaj- They satisfy 
(pphajiO) = —^j and remove the bad part of (pQ. 

(2) Correction from amplitude. The amplitude 6o which pairs (pQ via 
(1.13) has the expansion 60(1") >i ^o + X^'^li 'T^'''^ji where aj is a space 
function defined by ao (see. Proposition 4.2). Hence, the principal 
part of the amplitude of u^ [principal amplitude of u^) is 



^'o 



ffT 



l-t' l-t 



Go 



X 



l-t 



-^ T 






i=i 



t 



l-t 



(2.5) 



for eT/(l— t)^l. In particular, at t = we have 

fc-i 
60 (e^,x) ~ ao(x) + J^e"%i + ©(e"'^). 

i=i 
The principal amplitude converges to the given initial data for all 
a > 0, and so one might expect this is harmless. However, it is un- 
derstood that, when we try to get the pointwise estimate of solution 
via Grenier's method, we must take e^-term of the initial amplitude 
into account because it affects (implicitly) the approximate solution 
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at leading order (see, [2, 9, 18]). Therefore, we must remove e^-'flj 
for all J ^ 1 such that aj ^ 1 from (2.5), otherwise the approximate 
solution will differ outside a small neighborhood of t = 0. To do 
this, we construct &ampj as a solution to a kind of j-th linearized 
system of (1.11) with the condition 6ampj(0) = — Oj. At that time, 
there appear a phase correction (/>ampj associated with 6ampj via the 
system which 6ampj solves. 
(3) Correction from interaction. The third one comes from the structure 
of (1.11). As stated in introduction, the problem boils down to de- 
termining the asymptotic behavior of the solution (a , (p ) to (1.11) 
up to 0{h^). Suppose that the solution has two terms {h^^bi, /i^^^i) 
and {W^h2^ hP'^(p2) in its asymptotic expansion as h ^ 0. Then, the 
quadratic terms in (1.11) produce nontrivial hP^^^'^ -terras. For ex- 
ample, Vcp'^-Va^ has the terms /i^i+P^v^i . V62 and hP^^P^V(j)2-Vbi 
in its expansion. Again by (1.11), this implies that {dth ,dt(j) ) (and 
so {b ,(j) ) itself) also contains a hP^~^P^-term in its expansion. Re- 
peating this argument, we see that (6,0) has /i^-terms for all p 
given hy p = Ipi + mp2 with integers l,m ^ such that I + m ^ 1. 
So far, we have already obtained two kinds of correction terms, i?^pha,j 
and (/'amp J- Therefore, they interact each other to produce the third 
correction terms (pintj- 

• The supercritical cascade of phase shifts. We are now in a position 
to understand the cascade of phase shifts phenomena in the supercritical 
case. With the above three kinds of correction terms, we can describe the 
asymptotic behavior of the solution before the final layer. Recall that, in 
the linear WKB case a > 1, the principal nonlinear phase is given by only 
one phase function (/)o(r, y), and the notion of boundary layer comes from 
its expansion with respect to r around r = 0. In this case, not only 4>q but 
also all (j)pha.,i, 0amp,j5 4'mt,k produce a countable number of similar boundary 
layers. Thus, the cascade of phase shifts involves much more phase shifts 
and boundary layers than the linear WKB case. 

• One more correction at the final boundary layer. So far, the 
second difficulty listed in the beginning of this section is solved by three 
kinds of correction terms (j)pha.,i, <Aamp,j, and (pintj- We can describe with 
them the asymptotic behavior of the solution before the final layer. To give 
it also on the final layer, we need one more correction term, i;^equ) as in the 
case a > 1. This correction term defeats the first difficulty. This solves 
(1.18)-(1.19) and so it is independent of a. When a = 1, (/)cqu changes in 
to a solution of (1.18)-(1.20) (see Theorem 1.3). We next address why we 
need a modified (pcqu- 

• Resonance of correction terms and the nonlinear WKB case. 

Some of the correction terms associated with 4'pha.,i, (/"ampj, 4'mt,k, and (/)cqu 
may have the same order. This phenomena is a kind of resonance. The 
nonlinear case a = 1 is the simplest example. In this case, it happens that 
</'pha,i and (;^equ have the same order. Recall that if a = 1 then we work with 
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the modified (j)cqu solving (1.18)-(1.20), 



drbcqu + V0cqu " V6o + VcpQ ■ V6cqu + -bcqu^4'0 + -bo^Cpcqu = -^'^h 



dr4>cqu + V(^o • V(/)cqu + At'^ ^ {\y\ '^ * 2 Re 6o6equ) = 0, 
V^cqu|r=0 = Oi <^cqu|r=0 = K\y\~"^ * Wol'^) / il " 1) = "V'l 

(see Theorem 1.3). We will see later that 0pha,i is the solution to 

' 9^fepha,l + V0pha,l • V6o + V0O • V6pha,l + ^ Via,l^</'0 + 2^oA(/>pha,l = 0, 

^r-^pha,! + V(^o • V0pha,l + AT^~^(|y|"'^ * 2 Re 6o&pha,l) = 0> 
V&pha,l|r=0 = 0) 0pha,l|T=O = Klvl''^ * |aOp)/(7 " 1) = "V^l 

(see Remark 5.3). Let 0cqu be a solution to (1.18)-(1.19). One can check 
from above systems that (^cqu = 4'cqu + 0pha,i- Therefore, the modified cor- 
rection term c/^cqu is nothing but the superposition of the (usual) correction 
from equation 0cqu and the correction from phase (/'pha,!- 

3. Preliminary results 

3.1. Properties of the Y^ JR."-) space. We first collect some facts about 
the space Y^ defined in (1.16)-(1.17) (see, also [1, 9]). 

(1) i7,^cy;,2,ifsi^s2. 

(2) ^;,, = >;'[mi„(,,2-),oo] and so Y^^^^ C F^;,, if qi ^ gs- 

(3) If g < n then F/^^ = Y^tnin(p,g*),ooW I* implies l^,,, C Y^^^^ for 
Pi ^ P2 under q < n. In particular, Yp^j C Yr2** ^^i r2. ool ^^ ^ ^ ^' 
where 2** := (2*)* = 2n/{n - 4). 

(4) If n ^ 3 then any function / G X'' is written uniquely as / = 5 + c, 
where g £ Y^ 2!= ^2* 2) ^^'^ c is a constant. 

The first property is obvious by definition. The others follow from the 
following lemma which is a consequence of the Hardy-Littlewood-Sobolev 
inequality and found in [19, Th. 4.5.9] or [16, Lemma 7]: 

Lemma 3.1. // 99 G p'(IR") is such that V99 G LP(R") for p G]l,n[, then 
there exists a constant c such that (^ — c G L'?(]R"), with 1/p = 1/g + 1/n. 

We take a function / G Y^^^. Then, the indices p and q almost indicate 
the decay rates at spacial infinity of the function / and its first derivative 
V/, respectively. The second property means that V/ is always bounded 
and decays at the spacial infinity so fast that V/ G L^ . What the third 
property says is that / has a similar decay property. It can be said from 
the fourth property that 1^ 2 = ^^ in a sense, provided n ^ 3. Note that 
every g G Y^2 satisfies 5 ^ as |x| ^ 00 by definition (1.16)-(1.17). On 
the other hand, elements of X'^ do not necessarily tend to zero at the spacial 
infinity. We also note that if g = 2 then we have another definition of Y^2- 

(3.1) Yl^.i^'^) = LP(M") n X"(M"), 

which makes sense for s > n/2. 
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3.2. Basic existence and approximation results. Operate V to the 
equation for (j)^ in (1.11) and put w^ := \7(j)^. For our further apphcation, 
we generahze the system sHghtly. Let 

(3.2) Qi{a,v) = -{vV)a- -aV -v, 

(3.3) Q2{vi,V2) = -{vi ■ V)f2, 

(3.4) Q3(ai, as) = -\Vi\x\-^ * (aio^)), 
and consider a system of the fohowing form: 

' dtb^ = c5'Ql(6^u;^) + Ql(S^^i;^) + Ql(6^ W^) + R'^ + ir^A6^ 

(3.5) I dtw^ = 4Q2{w\w^) + Q2{Wlw'')+Q2{w'',W^) 

(3.6) 6f,=o = b^ wf,^, = wl 

where b takes complex value and w takes real value. Other notation will 
be made precise in Assumptions 3.2 and 3.3, below. In [9], the existence of 
a unique solution for this kind of system is shown with explicit coefficients. 
We shall summarize the parallel result. 

Assumption 3.2 (initial data). Letn ^ 3 and max(n/2 — 2, 0) < 7 ^ n — 2. 

We suppose the following conditions with some s > n/2 + 1; The initial 
amplitude 6g G i7''(]R") and the initial velocity Wq G ^n\ i^"^) for some 
qo G]n/(7 + 1), n[, uniformly for h £ [0, 1], that is, there exists a constant C 
independent of h such that Wb^WH" + ll^ollys+icron-i ^ C. 

Assumption 3.3 (coefficients). Let c^ and C2 be complex constants bounded 
uniformly in h, and let r^ be a real constant bounded uniformly in h. Suppose 
for some T* > and s > n/2 + 1 that f^ is a real-valued function of 
time and f^ G L^{{0,T*)); B^ and Ri are complex-valued functions of 
spacetime, and B^ G Li((0,r*); i7"+i) and R\ G L^{{Q,T*);H'); W^" and 
i?2 o,fe real-valued functions of spacetime, and W^ G L^((0, T*); y^"^2 ) and 
i?2 G -^^((0)^*))^''^ )■ Moreover, suppose all above functions are bounded 
in the corresponding norms uniformly with respect to h. 

Assumption 3.4 (existence of the limit). In addition to Assumptions 3.2 
and 3.3, we suppose the existence of limits of all b^, Wq, c^, r^, f^, B^, 
W^ , and R^ as h ^ in the corresponding strong topologies. These strong 
limits are denoted by bo, wq, Ci, r, f , Bi, Wi, and Ri, respectively. 

Proposition 3.5. Let Assumptions 3.2 and 3.3 be satisfied. Then, there 
exists T > independent of h, s, and qq, such that for all h G [0, 1] the 
system (3.5) -(3.6) has a unique solution 

{b\w'^)ec([0,T];H'xY^'+i 



Moreover, the norm of {b^,w^) is bounded uniformly for h G [0,1]. //, in 
addition. Assumption 3.4 is satisfied, then the pair {b^,w^) converges to 
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{b,w) := {b^,w^)\h=Q in C([0,T]; i^^-^ ^ y^-^) as h ^ 0. Furthermore, 
{b, w) solves 

( dtb = ciQi{b, w) + Qi{Bi,w) + Qi(6, Wi) + Ri + irAb, 

dtW = C2Q2{w, w) + Q2{W2,w) + Q2{w, W2) 

+ fit) {c2Q3{b, b) + Qs{B2, b) + Qsib, B2)) + R2, 

b\t=Q = bo, wit=o = Wo- 

Proof. The key is the following energy estimate for s > n/2 + 1 

+ {\\B'1\\hs+i + 11^2 Ik- + W'^W^Wh^ + \\VW^\\hs+i)E'' 

where E^ := Wb^Wjjs + ||Vw^|||^s. For more details, see the proof of Propo- 
sition 4.1 in [9]. D 

Remark 3.6. We intend to apply this proposition to the system (1.11). In 
that case, / (t) corresponds to f^"^, which is singular at t = if 7 < 2. 
Since f^{t) is only supposed to be integrable, we will see that the system 
(1.11) (and so the equation (1.8)) has a unique solution for 7 > 1 while the 
singularity. We also note that this corresponds to the fact that the Hartree 
nonlinearity is short range when 7 > 1. 

Remark 3.7. In the convergence part of Proposition 3.5, it can happen that 
n/2 + 1 ^ s — 1 > n/2. In this case, we use the definition (3.1) instead of 
(1.16). 

We conclude this section with a lemma which we use for the construction 
of a function (f)^ from the corresponding solution w^ = Vfj)^ to (3.5)-(3.6). 
This lemma is a consequence of Lemma 3.1. 

Lemma 3.8. // f satisfies |c^| ^ as |x| ^ cxd and Vf £ Y^2 f^f some 
s > n/2 and some q < n, q ^ 2* then tp G Y^, 



s+l 



1 ,<i ■ 



4. Proof of Theorem 1.3 

4.1. Strategy. In this section, we illustrate the strategy of the proof of 
Theorem 1.3 rather precisely. As in Section 1, we first introduce the semi- 
classical conformal transform: 

(1.7) u'{t,x) = -—^-—r(-^.-^\eicp(i '^' 



(l_t)n/2^ \l_t'l-tj'-\2e{t-l 
Putting r := ei /(I — t) and y := x/{l — t), we find that (1.3)-(1.4) becomes 
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Put h = e^^'^h and denote V'^ by ^'*. We note that e ^ is equivalent to 
/i ^ as long as a < 7. Thus, our problem is reduced to the limit /i ^ of 
the solution to 

(1.8) i/^a,V' + ^AV" = AT^-2(|y|-^*|^'^|2)V;^ # ^^(y) = ao(y). 

Our strategy is to seek a solution ^^ to (1.8) represented as 
(1.10) ^'*(r,2/) = a'^(r,|/)e*^''(^'^)/\ 

with a complex-valued space-time function a^ and a real-valued space-time 
function (/>''. Note that c^ is expected to be complex- valued, even if its initial 
value ao is real-valued. Substituting the form (1.10) into (1.8), we obtain 



a QA + i:\^r\ + ^r^'WvV^ * \a 



h\2- 



+ ih ( dro^ + (V/ • V)a^ + ^a^Ac/.'^ - i-^Aa^ 



0. 



To obtain a solution of the above equation (hence, of (1.8)), we choose to 
consider 



(1.11) 



dra^ + V^'^ • Va^ + -a^A0^ = i-Aa^, 



9,0'' + -IV^'^I^ + ArT-2(|yr^ * la'^p) 



0, 



=/i-7- 



ao, 0, 



=/i^- 



0. 



The point is that this system can be regarded as a symmetric hyperbolic 
system with semilinear perturbation. In Section 4.2, we first prove that 
it admits a unique solution with suitable regularity (see Proposition 4.1), 
hence providing a solution to (1.8) and (1.3)-(1.4). 

By (1.10), in order to obtain a leading order WKB type approximate 
solution it suffices to determine 0{h^) and 0{h^) terms of (f)^ in the limit 
/i — > 0. Letting /i = in (1.11), we formally obtain the 0{hP) term (6o,i^o) 
which solves 



(1.13) 



dthQ + V0 • V6o + 2&oA0o = 0, 

dt<Po + ^|V./.o|' + At^-2(|y|-^ * |5oP) = 0, 



^ Oo|r=0 — flO, <P0|r=0 







introduced in Section 1. The difficulty of finding /i^-terms lies in the follow- 
ing two respects; firstly, the equation (1.11) depends on h through the term 
i2Aa ; and secondly the initial data of (1.11) is moving at a speed h"i-°' . 
In Section 4.3, we give the time expansion of (60, 0o) around t = 0. We will 
obtain an expansion of the form 



(4.1) 6o(r,y) x J]r^%,(y), 0o(r,y) x J]r^J-V.-(y) 

j=0 j=l 

(Proposition 4.2). This expansion is essential in handling the moving initial- 
data. 
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In Section 4.4, we finally determine 0{h^) term in the case q ^ 1. What 
to show is the existence of the limits 

a'ir)-boir) ^\r) - Mr) ^ , .. 

7 ^ Ocqu(,TJ, > <Pcqu(TJ 

as /i ^ 0. A formal differentiation of (1-11) with respect to h suggests that 

(6equ5 0cqu) may solve the linearized system 

(1.18) 

1 1 i 

,9,0cqu + V</.o • V^cqu + ArT-^dx^^ * 2Re6;;5equ) = 0. 
By means of (4.1), the following estimates hold at the initial time: 



■h/.. , „ 

0{h'y-°' 



a^(hi-") — bQ{h'i-") boi^hi-") — ao "t i^ 



h h 

(j) (hi-") — (poihi-") (j)Q(hi-") aji-i) I 

— = 0{h ~-" ' 



h h 



■~f — a. 



Note that -^2_ > i for all a ^ 1, however, ^^ — '- > 1 if a > 1 and 




^2- > 1 for a" ^' ^ 1 >.""'"^r°->- "(t~^) ^ 1 if ^, ^ 1 OT^^ "(t~i) 

C— a 

if a = 1. In particular, 

(/)"(/lT-") — 0o(/l' 
h 

where (pi is defined in (4.1). Therefore, the OQi^) term is described by 
(6equ,</'equ) solving (1.18) with (1.19) if a > 1 and with (1.20) if a = 1. In 
Section 5, we consider the case a < 1. In this case, the above powers -^^ and 

^\^_^ are less than one, in general. Therefore, there appear several terms 

which is order less than 0{h^) in the expansion of (a , (j) ). We determine all 
these terms and obtain the asymptotic behavior of {a\(j) ) (see. Theorem 
5.1). 

In Sections 4 and 5, we mainly work with v = V0 instead of (jr itself. 
Note that, by means of Lemma 3.8, it is easy to construct (p from v . 

4.2. Existence of phase-amplitude form solution. According to the 
strategy in Section 4.1, we first show that the system (1.11) has a unique 
solution. 

Proposition 4.1. Let Assumption 1.2 he satisfied. Assume < a < 7. 
Then, there exists T > independent of h such that, for all h G (0, 1], there 

exists a unique solution ip G C{[hi-" ,T + hi-"]; H°°) to (1.8). Moreover, 
^h j^g written as 

where 

a'' G C([/i^,T + /i^];//°°)nC°°((/i^,T + /i^];i7°°) 
and 

<A'*GC([/i^,T + /.^];y(-/^^^]_(„/(^+^^^]) 



n c°^ ((/i^^ , T + /i^^] ; y.^ 



(n/7,oo],(n/(7+l),oo] / ' 
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Moreover, there exists a limit {bo,cl)o) := {a^,cl)^)\h=o belonging the same 
function space as (a'*, (j)^) (h > Oj, and (a'*, (f)^) converges strongly to {bo, cpo) 
as h ^ 0. Furthermore, (60, 0o) solves (1.13). 

Proof. We set velocity v^ = Vc/)^. Then, the pair {a^,v^) solves 

\ drv^ = Q2{v\v'')+T"'~^Q3{a\a''), 

(4.3) af <^ = ao, v'' <^ = 0, 

where Qi, Q2, and Q3 are defined by (3.2), (3.3), and (3.4), respectively. To 
fix the initial time, we employ the time translation r = t + /it-« , Then, the 
equation is 

idta^ = Qi{a\v^) + i^Aa\ 

(4.4) I 2 

I dtv'' = Q2{v\ v^) + (t + h—y-^Qs{a\ a'^), 

(4.5) af,=o = ao, vf,^, = 0. 

Now, the assumption 7 > 1 implies that (t + h^-'^)'^~ is integrable over 
(0,T*) for some T* > and its integral is uniformly bounded with respect 
to h. Fix s > n/2 + 1. Then, applying Proposition 3.5 with q = C2 = 1, 
r'' = h/2, f^{t) = (t + /i^)T-2, and Bf = W^ = E!} = 0, we obtain the 
existence time < T ^ T* independent of e and the unique solution 

(a^.^)GC([0,T];//^xy(;+(;^^^^,^,) 

to (4.4)-(4.5) such that ||a'*|||^s + llVf'^lll^s is bounded uniformly with re- 
spect to t G [0, T]. The upper bound depends only on HaoH/fs. Notice that 
Assumption 3.2 is satisfied for all s > n/2 + 1 and go £]^/(7 + 1); ^I- 

By the equation and the Hardy-Littlewood-Sobolev inequality, dtv'^ be- 
longs to C([0, T]; y(n/('Y+i) ooi 2)' Hence, we see from Lemma 3.8 that dt(j)^ G 
'"+^ ^ '2--" ^^ - 0, we also have 



^([0'^];^W;oo],(n/(,+i),oo])- Since </>f^o 
0'^GC([O,T];y.fi 



(n/7,co],(ra/{7+l),oo] '' 

Therefore, we have ^'^ ^ as |x| ^ 00. Then, again applying Lemma 
3.8 to v^, we conclude that 0^ G C{[0,T];Y^^^^ ^^ (n/{7+i) oo])- ^ince the 
existence time T is independent of s, we have a^ G C{[0,T];H^) and 
0^^ G C([0,r];y-/^^^] („/(^_^,)^^]). The bootstrap argument gives the C°° 
regularity with respect time. The existence of the limit (6o)</'o) and the 
convergence {a^,(f)^) -^ {bo,(j)o) as /i ^ follow from the latter part of 
Proposition 3.5 and Lemma 3.8. D 
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4.3. Time expansion of the limit solution near r = 0. By Proposition 
4.1, the system (4.2)-(4.3) has a unique solution even if /i = 0. We keep 
working with v^ = V(/)^ instead of (j)^ Write {bo,wo) := {a^,v^)\h=Q. The 
main difficulty of describing the asymptotic behavior of {a^,v^) as h ^ 

a 

comes from the fact that the initial data is given at r = hi-" . In order to 
handle this /i-dependence of the initial time, we give a time expansion of 
{bQ,wo) around r = 0. 
Note that (60,1^0) solves 

(4.6) drbo = Qi{bo,wo), drWo = Q2{wo,wo) + T"'~'^Q3{bo,bo), 

(4.7) 6o|^=o = ao, wo\r=o = 0, 
where the quadratic forms Qi are defined by (3.2)-(3.4). 

Proposition 4.2. Let {bQ,wo) = (bQ,'V(l>Q) be the unique solution to (4.6)- 

(4.7) defined by Proposition 4-i- Then, it holds that 

J 

(4.8) bo{T,y) = Y,r^'a,{y) + o{T^') mH^, 

j=0 

J 

(4.9) wo{r,y) = Y^r^^-^v,{y) + o{r^'~') ^n l(^/(,+i),oo],2 

i=i 
as T ^ for all J, where oq is the initial data for bo, Oj and Vj are defined 
by 

«j = ^ X^ Q 1(0^1,^^2) 

for j ^ 1, ui = (53(00, ao)/(7 - 1), and 
1 



7i- 1 



^ Q2{vk^,Vk2)+ X] Q3(Ofci,afc2) 

ki^l,k2^1,ki+k2=j ki^0,k2>0,ki+k2=j-l 

for j ^ 2 with the quadratic forms Qi defined by (3.2) -(3.4). Moreover, (J)q 
is expanded as 

J 
(1.24) Ur,v)=Y.^''~'^M + o{r^'-') ^nY^,^^^^^^^,^^_,,^^^^ 

i=i 
as r ^ for all J ^ 1, where ipj is given by cpi = j— ;(|y| '^ * l^op) o,i^d 



ipj 



i-7i 



Y^ ^(Vv^fei • V^JfeJ 



2 



ki^0,k2>0,ki+k2=j-l 

Remark 4.3. Using the "x" sign defined in Notation 1.1, the above three 
expansions (4.8), (4.9), and (1.24) can be written as 60(7") ^ '}2i=o'^^''^j ^^ 
H^, woir) X ZU^'"-'vj in l^(^/(,+i),^],2' ^^^ Mr) ^ EU^'"-'^j m 
^W7,oo],(n/(7+i),oo] ' respectively 
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Proof. We first note that, by the definitions of Qi, it follows for all s > n/2+1 
that 

\\Qi{b,v)\\Hs ^ Cs \\b\\Hs+i \\v\\y'>+-^ 

(n/(7+l),oo],2 

\\Q2ivi,V2)\\Ys ^CsWviWys \\v2\\ys+i , 

-'(n/(7 + l),oo],2 -'(n/(7 + l),oo],2 ^ (n/(7+l) ,00] ,2 

\\Q3ibi,b2)\\Ys ^ CsWhWHs IMIhs- 

(n/(7 + l),oo],2 

Therefore, we see that ai and vi are bounded in H°° and YP^,, -^-^ , ^, 
respectively. For simplicity, in this proof we denote Yf^,^ j^-^^-, , ^ by Y^. 

Denote J2j=o''''^''^J ^^'^ J2j=i'^''''^^'^j W 5; and vi, respectively. Then, it 
suffices to show that 

(4.10) 11^0 -SdlLoc([o,,].^oo) = o{t''^) yi ^ 0, 

(4.11) \\wo - ^dlL°c([o,r];y-') = o(r'^'"^) \/l ^ 1. 

step 1. Since bo G C([0,T]; iJ°°) and 6o(0) = ao = So, (4.10) is trivial if 
I = 0. We show (4.11) for / = 1. By the second equation of (4.6), it holds 
for s > n/2 + 1 that 



\\wo{T)\\ys^Ci \\wo{t)\\Ys dt + C2 t^~^dt, 
Jo Jo 

where Ci depends on s and C([0, T]; Y°°) norm of wq, and C2 depends on s 
and C{[0,T];H°^) norm of 60 • The right hand side is monotone increasing 
in time, hence this gives 

ll^o|lLoo((o,T-];y=) ^ C'lT ||wo|lLoo((o,i];ys) + C2t'^ ■ 

Choose r so small that Cit ^ 1/2. Then, we obtain 

lho|lLoc((o,^];yoo) =0(r'^"^) 
since s > n/2 + 1 is arbitrary. Again by the equation, it holds that 



wo-vi= Q2{wo,wo)dt+ €' {Q3{bo,bo- ao) + Q3{bo-ao,ao))dt. 
Jo Jo 

Since wq is order 0{t'^^^) in L°°((0,r]; y°°), the first integral of the right 
hand side is order 0{t'^'^^^) in L°^((0, r]; Y°^). Similarly, the fact that bo—ao 
is order o(l) in L°^((0, r]; i7°°) shows the second integral is order o(r'*'~^) 
in L°°((0,r];y°^), which proves (4.11) for / = 1. 

Step 2. We prove (4.10) and (4.11) for large / by induction. By the defini- 
tion of aj, an explicit calculation shows 

drbo = Qi{bo,wo) = Qi{bo,wo -n) + Qi{bo -Scr'^'^fi) + 3^(r'^ai) 

= Qi{bo,wo -V2) + Qiibo -ao,T'^'^'^V2) + Qiibo - ai,r'^"^fi) 

+ dr{T^ai+T^^a2) = --- 



i-i I I 



Qi{bo,wo-vi) + y]Qi{bo-ai,,T^^'-''^-\i_i,) + dr yr^'a, 



CASCADE OF PHASE SHIFTS FOR HARTREE EQUATION 23 

Similarly, it holds that 

drWo = Q2{wo,wo) + r'^'^Qsibo, bo) = ■■■ 

I 
= Q2iwo,wo-vi) + ^Qs (wo - Vi^,T^^^'^^^'^'^-'^Vi_i^+i 

h=l 

+ r'-^ Qz{boM-ai) + Yl Qz{bo-ai,y^'~'^\i-i,) 
\ h=o 

(i+i 

Integrating these identities with respect to time, we obtain 

(4.12) bo-ai = J {Qi{bo,wo-vi) + Y.Q^(''o-^h,t''^'~''^~\-h)\dt 



and 
(4.13) 



I 



Wo - vi+i = I I Q2{wo,wo -vi) + Yj Q-i (^0 - vi^,f<^^ 'i+^) ^Vl-l^+ij dt 

+ f f'~^ ( Qz{boM-ai) + Y, Q3{bo-ai„t"'^'-''^ai_i,) j dt. 

Now, let L^ 1 be an integer. If (4.10) holds for / ^ L - 1 and (4.11) holds 
for / ^ L, then we see that (4.12) gives (4.10) with I = L. On the other 
hand, if both (4.10) and (4.11) hold for I ^ L, then we obtain (4.11) with 
l = L + l from (4.13). 

The expansion of (po is an immediate consequence of the expansion of 
Wo = V0O. Since 

Q2iVki,Vk2) + Q2iVk2,Vki) = -^{Vk^ -Vki) = -V f -Wfci -Vk^ + -Vk2 -Vkij , 



we deduce from the definition Vj that 



V^j = V 



J — ""J 



i-7i 



:V 



E 



:(V(^fci • Vv^fej) 



ki^l,k2^1,ki+k2=j 



ki^0,k2^0,ki+k2=j-l 



By Lemma 3.8, (pj belongs to YP^ 



(n/7,oo],(n/(7+l),oo] ' 



D 
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4.4. Asymptotic behavior of the phase-amphtude form solution. 

The following proposition completes the proof of the theorem. 

Proposition 4.4. Let Assumption 1.2 satisfied and q ^ 1. Let T > and 

{a^,v^) be as in Proposition 4-1- Let (6o,^fo) := {of^ ,v^)\h=Q- Then, there 
exists (6equ,ifequ) G C([0, T] ; i7°° X ^('^/(^+i)^oo],2) '^'"^^ ^^"''^ ^^^ following 
asymptotics holds: 
(4.14) 

a'^ir) = 6o(r) + /i6cqu(r - h^)+o{h) in C{[h^ ,T];H^), 
v\T)=w,{T) + hw,^,{T-hi^) + o{h) «nC([/i^,T];y(-/(^+,)^^]^2)- 

Moreover, (ftequj^i'equ) solves 
(4.15) 

drhcqu = Ql{bo,Wcqu) + Qlibequ,Wo) + i-Abo, 
drWcqu = Q2{wo,Wcqu) + Q2(li'cqu, ^^o) + T^~^ (Qsl&O, ^cqu) + Qsibcqu, ^o)) 

with the data 

\ if a > 1, 

(4.16) Ocqu|r=0 = 0, ^^^equ|r=0 = < •. i 

l^-fi ^/a = l, 
where vi is defined in Proposition J^.2. 
Remark 4.5. Since 

a a 

bcqn{r - h-'-" ) = bcqnir) + o(l), ti;cqu('r - h-y-^: ) = Wcqn{r) + 0(1) 
by continuity, (4.14) implies 

a'^ir) = 6o(r) + /ifeequ(T) + o{h) in C{[h^ ,T];H^), 
v\T)=wo{T) + hw,q^{T) + o{h) in C {[h^ ,T];Y^'^/^.^^,)^^^^,) . 

From this asymptotics and the transforms (1.7) and (1.10), we immediately 
obtain the asymptotics (1.21). 

Proof. Let {a^, v^) be the solution to (4.2)-(4.3). Let {bo,wo) be the solution 
to (4.6)-(4.7). We put 

. a'*(r + /i^=^,y) -6o(r + /i^^,y) 
b (r, y) := 



w^iT,y) :-- 



h 

V {t + /it-" , y) — wo{t + /iT-Q , y) 



h 

Then, {b^,w^) solves 
(4.17) 

d,t = hQi{b^, w^) + Qi(6o, w^) + Ql(6^ wo) + i^A6o + ^^A6^ 
S^w;'^ = /^Q2(^^ w^) + Q2(w'o, ti''') + Q2(w'\ ti'o) 

+ (^ + /i^)7-2 (/jQ3(6'*, ft'*) + Q3(6o, ft'^) + Q3(6\ 60)) , 

, . h _ bo\r=0-Hh'^) h _ Wolr=O^Woihf^ 
l4.i»j 0|^^o - , U)\r=0 - -f^ ' 
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We apply Proposition 3.5 with these initial data and Ci = C2 = h, r^ = h/2, 
fh(t) = {t + /i^)T-2, B'^ = ^2 = ^0, W^ = W^ = wo, R\ = {i/2)Abo, 
and i?2 = 0- Note that the initial data (4.18) is uniformly bounded if a ^ 1 
since an application of (4.8) and (4.9) gives 



h 



0{h f-Q 7-a) in H"^ 



woi^=o-woihf—) 7(^-1) . 

— U{n T-a ) m Yf 



00 



The term R\ satisfies ||-Ri Hh^ ^ ||6o||h»+2/2- Therefore, if s - 2 > n/2 + 1, 
that is, if s > n/2 + 3 then Proposition 3.5 provides the unique solution 

{h\w^) G C([0,T - hi^Y^H^-^ X y(;7(\+i),^],2) for h G [0,1]. More- 

over, {b^,w^) converges to {b,w) := {b^,w^)\h=Q in C([0,T — h^^^"]; H^~^ x 

'^Ui/( +1)001 2)' ^^ follows from (4.9) that lim/j^o^u-=o = if ct > 1 a-iid 

lim/i_^o ii'|^=o ~ ~^i if a = 1. Hence, {bequ,Wequ) solves (4.15)~(4.16). D 



5. Supercritical caustic and Supercritical WKB case 

5.1. Result. In this section, we treat the case q < 1 < 7. As presented in 
Section 4.1, the asymptotic behavior of the solution (1.3)~(1.4) boils down 
to the asymptotic behavior of the solution to (1.11). By means of Lemma 
3.8, we work with (a^,^'') := (a'', V(/)^) which solves (4.2)-(4.3). 

The main difficulty lies in the fact that the initial data (4.3) is moving at 

a 

the speed /it-" (see Section 2.2). From the expansion (4.1) of (60,1^0) •= 
{a ,v )\h=Q, we deduce that {a ,v ) contains the terms of order 

a(7i — 1) ja-^ 

0{h 1-" ) and 0{h-y-") 

for all i,j ^ 0. Note that some of these orders are less than one if q < 1. 
This is the feature of the supercritical WKB case, and the problem comes 
from this point. Moreover, the above terms interact each other and there 

a^-yi—l) 

appear all the terms whose order is given by the finite combination of h t-^ 

and h~<-°' . Thus, we see that [a ,v ) contains all the terms whose order is 
written as 

0{h f— ), 0^/2^^1. 

For our purpose, we determine all these terms up to 0{h^). Therefore, it is 
natural to introduce a set P defined by 

(5.1) P := I "(^^^-^^^ 0^/2^/1, ^ "^^^^ " ^^^ < 1 

[ 7 — a 7 — a 

Set N := jijP — 1, and number the elements of P as = po < Pi < • • • < 
Pn < 1. For any Pi^,Pi2 G P, either pi^ +Pi2 G P or pi^ ^pi^ '^ 1 holds. For 
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example, if 7 = \/3 and a = \/3/4 then, po = 0, 



VS-l a(7-l) 2(^/3-1) q(27-2) 



Pi = — ^ — = , P2 



3 7 — 0; 3 7 — Q 

, . V3 «7 /:j 1 a(37-3) 

(5.2j p3 = = , p4 = V3-l 



P5 



3 7 — a 7 — Q 

2^/3-1 a(27-l) 



3 7 — a 

and N = 5; and if 7 = 2 and a = 1/3, then po = 0, 



1 ah - 1) 2 a(27 - 2) 07 

Pi = 7 = ' P2- - - 



(5.3) p3 



5 7 — a 5 7 — Q 7 — a 

3 _ a(37 - 3) _ 0(27-1) 
5 7 — a 7 — a ' 



4 a(47 - 4) a(37 - 2) 2a-f 

P4 = - ~ ~ ~ 



5 7 — a 7 — a 7 — a 

and A^ = 4. 

To state the result, we also introduce several systems. Let Qi, Q2, and 
Qs be quadratic forms defined in (3.2), (3.3), and (3.4), respectively. Let ai 
and vi be sequences given in Proposition 4.2. Then, for any ^ i ^ N, we 
introduce 

drbi= ^ Qi{bp^,Wp^), 

V3+Pk=Pi 

drWi= ^ {Q2{Wp^,Wp^})+T'^''^Q^{bp^,hp^)) , 

Vj+Pk=Pi 



(5.4) 



(5.5) 



, / N I —11 if 3/ such that pi = -z-^, 
6i(0) = < n ^-a^ 

I otherwise, 

\-vv if 3/' such that pi = °'Z^'~^ , 
Wi[U) = < ^ " 

otherwise. 



We also introduce a system for (6cqu,ii'cqu) 
(5.6) 

<9r6equ = Ql{bo,Wcqu) + Ql (&cqu, it^o) + ^ Ql{bp^,WpJ + 2^^0, 

Pi+Pfc = l 

C^rti'equ = Q2{wo,Wcqu) + Q2(li'cqu, tt^o) + T^"^(<53(^0, &cqu) + Q3(^cqu, ^o)) 

+ Yl {Q2iWp^,Wp^)+r^~'^Q3ibp,,bp^)) , 
Pj+Pk=i 



—ai if 3/ such that 1 = -22_ 

t ^ ,-v ' 

otherwise. 



6cqu(0) = ^ "' -^-^ 



(5-7) 

J-t-r if3/'suchthatl = ^2L^, 

^cqu(O) = < . ^ " 

(J otherwise, 
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where (60, ^^o) is the solution of (1.13). If there is no pair {j,k) such that 
Pj +PA; = 1, we let X]pj+pfe=i = 0- It ™ay happen (see (5.2)). 

Theorem 5.1. Let assumption 1.2 he satisfied. Assume < a < 1. 
Let P he as in (5.1) and N = (JP — 1. Then, there exists an existence 
time T > independent of e. There also exist {bj,4>j) £ C{[0,T];H°° x 
^K7,oo],(n/(7+i),oo]) CO ^ i ^ N) such that {bi,Wi) := {bi,V(t>i) solves 
(5.4)-(5'.5), and (6equ, c^cqu) G Ci[0,T];H°- x l^(^/^,„o],(n/(7+i),oo]) «^c/^ that 
{bcqu,Wequ) '■= (^cqu, V</>cqu) solvcs (5.6) -(5.7). Morcover, the followings 
hold: 

(1) (/.o(t) X Ej°li^^-'^Vj (in the sense of (1.24);. 

(2) The solution u^ to (1.3) -(1.4) satisfies the following asymptotics for 
all s ^ 0; 

(5.8) 



sup 

tG[0,l-T-ie«/T] 



r\' u%t)e 



-i<I>=(i) 



(1 - t)«/2 



^£(i)e*2iiF^ 







L2 



as £ ^ with 



{i-f)p. 




tCT X 



1-t' 1-t 



(5.10) A''(t,x) = 6o 



ET 



1-t' 1-t / ^ \ ""^" V^-* i~ v/ 

Remark 5.2. In (5.9), the time variable of (pj {j ^ 1) is not £-' /(I — t) but 



feT 



ST 



1-t l-t 



ET . 



Although this variable is not stable on the final layer 1 — t = T~^e^, this 
choice is suitable when we work with the well-prepared data (see Section 
5.2). Of course, the Taylor expansion 



teT 
l-t 



Y.{-e^f{d1 



fc=0 



l-t 



will exclude the variable teT/(l — t) from $^, however, we do not pursue 
this point any more. 

Remark 5.3. Let us classify the phase functions in (5.8) according to the 
notion in Section 2.2. If 0i(O) ^ (resp. 6j(0) ^ 0), that is, if there 
exists a number / ^ 1 such that pi = ""''_"" (resp. pi = -^5~)) then 0j is 
the correction from phase (resp. correction from amplitude); in particular 
4>i = 4>pha,i (resp. (pi = (/>amp,/)- On the other hand, if (/>j(0) = 6i(0) = then 
(pi is the correction from interaction; in particular (pi = (pint,l' for some /'. 
Notice that the summation in the system (5.4) is decomposed as 



E 



E + E 



Pj+Pk=Pt {j,k)={i,0),{0,i) pj+p^,=pi,jk^O 
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The second sum is the interaction term, which is an external force. When 
4>i is the correction from interaction, it always has a nonzero interaction 
term. Otherwise, (pi = since the system for (j)i is posed with the zero 
initial condition. There is a possibility that the correction from phase (resp. 
correction from amplitude) has an interaction term Indeed, it happens if 
there is a triplet j,k,l such that pi = ^^^ = Pj +Pk (resp. pi = ^^ = pj + 
Pk) and jk ^ 0. In this case, there is a resonance between the correction from 
phase (resp. correction from amplitude) and the correction from interaction. 
4>eqn solving (5.6)~(5.7) is the correction from equation with or without 
resonance: If i;^(0) ^ (resp. 6(0) ^ ) then there is a resonance with the 
correction from phase (resp. correction from amplitude); if Xlp +» ^ ^ then 
there is a resonance with the correction from interaction. We note that the 
resonance among the correction from phase, the correction from interaction, 
and the correction from equation (resp. the correction from amplitude, the 
correction from interaction, and the correction from equation) may happen, 
and that, however, the resonance between the correction from phase and the 
correction from amplitude never happens because there is no pair I, I' such 
that 7? = 7?' - 1 if 7 > 1. 

Proof. As presented in Section 4.1, the asymptotic behavior of the solution 
(1.3)-(1.4) boils down to the asymptotic behavior of the solution to (1.11). 
By means of Lemma 3.8, we work with [a ,v ) := {a jVcj) ) which solves 
(4.2)-(4.3). 

The existence of {a ,v ) and the expansion of (pQ are already proven in 
Propositions 4.1 and 4.2, respectively. Let P be as in (5.1), N = 'i^P — 1, and 
Pi G P {i = 0,1, ... , N) be such that {pij^o ~ ^ ^^^ Pi < Pi+i- It suffices 
to show that {a^,v^) is expanded as 

N 



a^{T + hi-" ) = bo{T + hi-'-) + Y^ hP'biir) + /i6cqu(T) + o(/i), 

N 

v^{t + /i^^^) = wq{t + /i^^^) + ^ h^'"Wi{T) + hwequir) + o{h) . 



1=1 
Plugging this and 

a. a 

bcqnir) = 6cqu(T + h^-" ) + o(l), Wequir) = Wequ(T + h^-" ) + o(l). 

to (1.10) and (1.7), we obtain (5.8). 

Step 1. We first prove by induction that 

k 

a^ir + hi^) = bo{T + /i^) + Y, h^'biir) + oih^"), 
(5.11) 



'j'^{t + h-i-") = wq{t + h-i-") + Y^f^^'Wiir) + o{U 



4 = 1 



holds for k = N, where (hi, Wi) is a solution of (5.4)-(5.5). One verifies from 
Proposition 4.1 that (5.11) holds if A; = 0. We put K G [1, A^]. We assume 
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for induction that (5.11) holds for k = K — 1, and put 



K-l 



b'}^{T) = h'P^ a^{T + hi^) - 6o(r + hi^) - Y, h^'biir) , 
By the equation for {a^,v^), we see that {b^,w^) solves 



d^w'k = hP^Q2{w'k,w1,)+Q2{W^,w1,)+Q2iw'kM) 

+ (r + h^y-^ (hP^Qsib'k, b\) + Q^{B\, b\) + Qj,{b\, B\) 



+ Rl 



K-l 



5^(0) = h-P^ ao - boih—) - Y, h''%iO) 

^^(0) = /l-P^ (-wo{h^) - Y hP^WiiO)] : 



j=l 



where 



K-l 



K-l 



i=0 



Ri= Yl Qi(.bn,WpJ 

Pi <K,Pk<K,pi +pk =pk 
Pi <K,Pk<K,pi +pk >PK 

ih^'PK ^ 

+ ^ Y ^'-^^^ 

1=0 

-^ Y Qi{bp„WpJ, 

Pi <K,pk<K,pi +pk =Pk 
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and 




Rl^- 


Pi<K,Pk<K,pi+pk =PK 



Pi <K,Pk <K,pi +pfc >pk 
Pi <K,Pk <K,pi +Pk >PK 
Pi<K,Pk<K,pi+pk =PK 

Moreover, applying the time expansion of (60, wq), we deduce that {bj^{0),Wj^{0)) 
is uniformly bounded and that, as h ^ 0, 

( ... ^ 

^ \-aj> if 3/ such that px = p^, 
1 otherwise. 



tt;^(0) = - h-P^ Wo{h^) - Y ^ 



a(7J-l) 

-1-°' Vj I , 



{r,^'^^<PK} 



_^ \-Vj>' if 3j" such that px = ° ?jLa , 
1 otherwise. 

Note that either 6^(0) ^ or w^{0) -^ holds for all K since 7 > 1 implies 
there is no pair (j, j') such that 7j = jj'—l. Therefore, we apply Proposition 
3.5 to obtain the solution {b^,w^). Put {bx^WK) '•= {b%-,w^)\h=Q. Then, 
it solves (5.4)-(5.5), and so (5.11) holds for k = K. By induction (5.11) 
holds for k = N. 

Step 2. Mimicking the argument in Step 1, we construct (6equj li'equ) such 
that 

AT 

a^{T + /i^) = 6o(r + /i^) + Y h^'biCr) + hbecinir) + o{h), 



i=l 

N 

a. a 

v"[ — ; _; ■ ' : 

1=1 



'j''{T + h-y-") = woir + h-r-") + ^ /iP»ti;i(r) + hwcquir) + o{h), 
holds. Notice that (^cqm^i'cqu) solves 



drbeqn = Ql{bo,Wcqu) + Qlibcqn,Wo) + ^ Ql{bp^,Wp^) + ^^^0, 
drWequ = Q2{wo,Wequ) + Q2('Wcqu, Wq) + T^~^(Q3(&0, ^u) + Qsibcqn, ^o)) 



Pj+Pfc=l 
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6(0) = < 7-« ' ^^(o) = <^ ' 7-" ' 

lO otherwise, lO otherwise. 

D 

Remark 5.4. By a similar proof, we obtain higher order approximation. Let 
< a < 7 and 7 > 1. We modify the set P defined by (5.1) as 

[7 — a 7 — a J 

and number the elements of P' as = po < Pi < • • • < Pk < ■ ■ ■ ■ Then, we 
have, for all k, 

k 
a^{T + /i^) = 6o(t + h^) + Y, h^'biir) + oih^"), 

(5.12) 

(p'^ir + /i^) = Mr + hi^) + ^ hP^cpiir) + o(/iP'=). 

i=l 

Plugging this to (1.10) and (1.7), we obtain higher order approximation of 
the original solution. It is important to note that the (5.12) has the same 
form in the both a ^ 1 and a < 1 case. When we concerned with higher 
order WKB-type approximation, we need four kinds of correction terms 
even if a ^ 1. From this respect, the supercritical WKB case a < 1 can be 
characterized as the the special case pi < 1. 

5.2. Well-prepared data and general data. We conclude this section 
with some remarks about the well-prepared data. By the semiclassical con- 
formal transform (1.7) and Grenier's transform (1.10), the leading order 
WKB-type approximation of the original solution u^ to (1.3)-(1.4) is re- 
duced to the approximation of the solution {a^, (p^) to 



(5.13) 



with 






(5.14) a!" _^ = ao, 0^ a = 0. 

up to order 0{h^). Note that (5.13)-(5.14) and (1.11) are the same. As 
shown in Proposition 4.1, there exists a limit (6o)0o) '■= {o.^,4'^)\h=o which 
solves (1.13). Now, we consider the distances 

d'iit) ■.=a\t) - bo{t), d%t) ■.=<t^\t) - Mt). 

If these distances are order o{h^) then we immediately obtain the WKB- 
type approximation bo exp(ieT cpo) of n^ (recall that h = e t ). However, 
unfortunately, the following two respects prevent us: The first one is the h- 
dependence of the equation (5.13), and the second one is the /i-dependence 
of the initial time (5.14). The first problem is handled by employing the 
correction term (6cqu, i;^cqu) solving (1.18) and, therefore we discuss about 
the initial data in the followings. 
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The given initial data (5.14) is written as 

dllih"^) :=ao - boQi'^), d'l{h'^) := - 0o(/i^). 

The main difficulty in the case q ^ 1 is the fact that these terms become 
larger than 0{h^) as /i ^ 0. The simplest way to overcome this difficulty is 
to modify the initial data (5.14) into 

(5.15) a!" .^ =bo(h'^), (j)^ ^ = 6n(h^). 

\T=h-i-°' \T=hi-" 

which ensures d^{h'i-°') = d^{h~i-°') = 0. Note that (5.13) with (5.15) is the 
same as (1.14). Back to the transform (1.7), this initial data corresponds to 
the well-prepared data 

I |2 
a ■ \x\ <y_ -1 a_ 

(1.6) tif^^Q(x) = 6o(e^;a;)e~*^^ exp(JeT (j)Q{e't ,x)). 

This initial condition is rather natural in the supercritical case a < 1, and 
the original initial condition in (1.4) is a kind of constraint (see, Section 2.2). 
If we use this well-prepared data, we do not have to consider any correction 
term other than (6cqu; 4'c(\u) and, for all < a < 7, it holds that 

a^{T) =bo{T) + /l6cqu(T) + 0{h), 

4>^{r) =Mr) + /i0equ(r) + o{h) 

with (6o,i;^o) and (6cqu5<Acqu) solving (1.13) and (1.18)-(1.19), respectively, 
and so that the asymptotic behavior of the solution u^{t,x) to (1.3)-(1.6) is 
given by 

This approximate solution is the same one as in the case a > 1. We still need 
0(/i^)-correction term (6equ) '/'equ) because it comes from the /i^-dependence 
of the equation which (a'*, </>'*) solves. 

In Theorems 1.3 and 5.1, we took another way. By the expansion of 
(60, 0o) around r = 0, there exist nonnegative integers ki,k2 depending on 
a and 7 such that 

fci 
daih'i-") =ao — h^Qi'i-") = — y^h-'~°'aj + o{h ), 

i=i 

d^{h^-'^) = — (j)o{h-i-°') = — y^h T-" ipj + o{h). 

We subtract the main part — ^jLi h''-°'aj and — Xliii ^ ''"" Vj by con- 
structing appropriate correction terms (correction from amplitude and cor- 
rection from phase, respectively). Indeed, if we let the correction terms 
{hj,4>j) (0 ^ J ^ N) and (6cqu, ';^equ) be defined as in Theorem 5.1 then, at 
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the initial time r = /it-" , it holds that 

N 

a\T) - boir) - Y, h^'^b.ir - h^) - /i6equ(r) 



i=i 



jr=/i7- 



/ 



bo{hi-'^' 



V 



^ h-f-'- aA+h (6equ(0) - bcqnih-'-" ) 



'€{'^0;^^!} 



N 



(f^'^ir) - Mr) - Yl h'^'Mr -h — )- /l,/.equ(t) 

1~^ I \n hl-a 



\T=h' 



(t>o{h~- 



a(-fl-l) 






o{h'] 



It is important to note that the time variable of (6j, (j)j) (1 ^ j ^ N) is not 

T but r — h'-'-" . This choice is the key for the above cancellation. By the 
transform (1.7), the time variable of (pj (1 ^ j ^ N) in the definition (5.9) 
of ^^{t) should be given by 



£ T « tet 



l-t ' l-t 
These correction terms allow us to work with general data. 

6. Proofs of Theorems 1.5 and 1.6 
Recall that the system we consider is 

( dripr''-^) + drilpvr'''^) = 0, 



(1.26) 



drV + vdrV - At" ^drVp = 0, 



n— 1 



drir'^-'drVp) = pr 

, P\r=0 = |aOp, l'|r=0 = 0, 



where r = \x\. We introduce the "mass" m and the "mean mass" M 



m{T, r) := M{t, r^ := r"-^a^$(T, r) = / p{t, s)s''''ds. 

Jo 

We also set r7io(r) := ?7i(0,r) and Mo(r) := M(0, r). Combining the first 

and the third equations of (1.26), we obtain 

(6.1) drin + vdrin = 0, 

where we have used {pvr'^~^)\r=o = 0. To solve this equation we also intro- 
duce the characteristic curve X{t, R): 

^ = v{r,X{T,R)), X{0,R)=R. 
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Denoting differentiation along this characteristic curve by / := d/dr^ the 
mass equation (6.1) and the second equation of (1.26) yield 

(6.2) m' = 0, 

(6-3) V = Ar"-^, 

(6.4) X' = V. 

We solve this system with the initial data 

(X,m,f)|^=o = (-R,"^o(-R),0), 

where -R > parameterizes the initial location. 

By (6.2), the mass m remains constant along the characteristics, that is, 
m{T,X{T, R)) = mo{R). Therefore, (6.3) and (6.4) yield 

AmoCR)r"-4 
This equation is studied also in [29]. 



(6.5) X" = -"'"'^Zi ' XiO,R)=R, X\0,R)=0. 



Proof of Theorem 1.5. It suffices to show that X{t,R) ^ holds for some 
R > and t ^T*. This argument is similar to that for the blow-up for the 
nonlinear Schrodinger equation by Glassey [17]. Since A < 0, we see from 
(6.5) that X" ^ 0, and so that X' ^ X'{0) = and X ^ X{0) = R for all 
r ^ 0. Therefore, again by (6.5), we verify that 



|A |moT" -^ 



X'\t,R) ^ - '^'3, = -\X\RMoiR)T 



n—4 



Note that Mq{R) > for some R provided oq is not identically zero. Now, 
fix such R. Integrating twice with respect to time, we obtain 

[n — 2)[n — 3) 

which yields X(t,R) ^ for large time. The critical time is not greater 
than 

(n-2)(n-3) \^ 

|A|sup^^o^o(i?)y 
since R is arbitrary. D 

Proof of Theorem 1.6. In order to clarify the necessary and sufficient condi- 
tion, we repeat the proof in [12]. We first note that dX/dR{0, R) = 1, and 
that the solution is global if and only if dX/dR{t, R) > for all r ^ and 
R^O. 

By (6.5), we have 

X" 



XniQ 



X3 ■ 
We multiply this by X' and integrate in time to obtain 

(Y'\2 n2 , -^"^0 Amo Amo „ 

where we have used (6.5) again. This yields 

{Xy = 2{X'f + 2XX" = ^^. 

R 
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Then, integration twice gives 
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X' 



R^ + 



Xm, 



,2 



i?2 



t\ 



Since X ^ R > from (6.5) and ^'(0) = 0, we see 
(6.6) 



X(r,i?) = Wi?2 + 



XniQ 



An explicit calculation shows that 



d^i o^ l + A(i|aop-Mo)r2 



VI + AMoT^ 

where mo = MqR'^ = J^ joops^ds by definition. Note that liuiji^Q Mq = 
|ao(0)p/4 since oq is continuous. Hence, dX/dR{t, R) > holds for all r ^ 
and i? > if and only if ^|aoP - Mq ^ for ah i? > 0, that is, 

(6.7) \ao{R)\^^—J \ao{s)\h^ds 

for all R > 0. Moreover, the critical time is given by 

1 
2 



XmaXr>o (2Mo(r) - |ao(r)|2; 
The condition (6.7) can be written as 

'mo(-R)" 



5i? 



i?2 



^0. 



We take i?o so that mo(i?o) > 0. Then, it holds only if 

moiR) > ^^R' 
Rq 

for all R ^ i?0) which fails if lim^^oo t^o{R) < oo, that is, if oq £ -^^ 



D 



Remark 6.1. If n = 4 then (1.26) becomes an autonomous system. Now we 
consider the classical solution of autonomous model 

' dripr"""^) + dripvr''-^) = 0, 

drV + vdrV — XdrVp = 0, 

drir^'-'drVp) = pr'''\ 

_ P\r=0 = PO, V\^=0 = VO, 

where {T,r) G R+ x IR+. It is shown in [25] that, under the assumption that 
n ^ 3, po £ L^((0, oo),r"^^(ir), vq{0) = 0, and Vo{r) ^ as r ^ cx), the 
corresponding solution is global if and only if A < and 



Voir) 



2|A| 
(n — 2)r' 



n-2 



Po{s)s'^ ^ds. 
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Remark 6.2. The function T = dX/dR is called the indicator function. As 
in the above proof, the solution blows up if and only if F takes non-positive 
value. Moreover, the solution is given by 

p{t,Xit,R))- ^°(^)^"" 



X"-i(t,i?)r(t,i?)' 

v{t,X{t,R))=^{t,R). 

Example 1.7 is easily checked by this form since the characteristic curve X 
is given explicitly by (6.6) in the case of A < and A^ = 4. 
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